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Abstract. The squeezing dynamics of a damped harmonic oscillator are studied 
for different types of environment without making the Markovian approximation. 
The squeezing dynamics of a coherent state depend on the reservoir spectrum 
in a unique way that can, in the weak coupling approximation, be analyzed 
analytically. Comparison of squeezing dynamics for Ohmic, sub-Ohmic and supcr- 
Ohmic environments is done showing a clear connection between the squeezing-non- 
squeczing oscillations and reservoir structure. Understanding the effects occurring due 
to structured reservoirs is important both from a purely theoretical point of view and 
in connection with evolving experimental techniques and future quantum computing 
applications. 
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1. Introduction 

It is known from the theory of open quantum systems that the environment affects 
the state of the system by inducing decoherence and heating, which very quickly 
destroy quantum superpositions and entanglement jlj. This effect is a major obstacle in 
quantum information science which relies on the fragile quantum properties. Very often 
it is assumed that the environment does not have a memory, i.e., changes in the reservoir 
induced by the system do not have a back action effect on the system. This assumption 
is known as the Markovian approximation. However it is not valid in general for strong 
couplings or for certain reservoirs with long memory. These include, e.g., photonic band 
gap materials [2] and atom lasers [3]. The need for non-Markovian theory depends on 
the properties of the environment. Recently, as experimental setups approach the limit 
where the Markovian approximation ceases to be valid, proposals for non-Markovian 
quantum computation have appeared in the literature [I]. 

When the non-Markovian effects are taken into account the evolution of the system 
is different from the Markovian case. To study non-Markovian effects I use the quantum 
Brownian motion model, describing a particle in a harmonic potential coupled to a 
quantized bosonic thermal reservoir [TJ. An exact solution to the master equation 
describing the time evolution of the system exists E] . Of course the analytical form 
of the solution depends on the spectral distribution of the reservoir. In this paper I 
consider three different types of spectral distributions, Ohmic, super-Ohmic and sub- 
Ohmic. This allows to compare the effects of different reservoirs on a given quantum 
system which correspond to different physical realization of, e.g., a qubit. In this way 
we can clarify the microscopic processes underlying the dynamics of exemplary open 
quantum systems. 

In this paper I consider an initially squeezed coherent state. The Heisenberg 
uncertainty principle for the variances of the dimensionless quadratures x and y states 
that (Ax) 2 (Ay) 2 > 1/4. Squeezed states are the ones where the variance of one 
quadrature is smaller than that of the vacuum, i.e., 1/4. The other quadrature will have 
a larger variance in order not to violate the uncertainty principle. Previously we have 
studied the effects of an Ohmic, super-Ohmic and sub-Ohmic reservoir to the dissipation 
of a quantum harmonic oscillator [7] and the decoherence of a coherent superposition 
state [8]. In this paper I investigate how the short-time non-Markovian dynamics of a 
squeezed state evolution is affected by these different reservoirs. 

Squeezed states of light can be utilized in quantum information processing [9j [10] 
and quantum communication [TT] to improve error rates. They have been used to 
construct entangled states and to demonstrate quantum teleportation of continuous 
variable quantum states [121 [13]. Squeezed states are also used in quantum metrology 
and high precision measurements [T4"l H~5] . To the best of the author's knowledge, 
squeezed states of the harmonic oscillator have been studied in the non-Markovian 
regime only for the Ohmic reservoir [16] , where it was shown that the squeezing dynamics 
are a result of virtual processes between the system and the bath. 
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The paper is organized in the following way. In section II we present the model 
and master equation for the system. Section III introduces the different reservoir types 
used in the comparison. The main results of this paper are given in section IV. Finally 
section V contains the conclusion. 



2. The system and the master equation 

Our model is the quantum Brownian particle in harmonic potential which consists of a 
quantum harmonic oscillator linearly coupled to a bath of quantum harmonic oscillators. 
The Hamiltonian of the system and the bath, in units of H, are Hs = u>o(a^a + 1/2) and 
He = ^2 n ^n{bn°n + 1/2), respectively, and the microscopic interaction Hamiltonian is 

Hl= 9 {a + a ^k n (b n + bi) } (1) 

* n 

where a(b) and a^(b^) are the annihilation and creation operators of the system (bath), 
g is a dimensionless coupling constant, k n gives the coupling between the system and 
each individual environment oscillator, and Uq and u n are the frequencies of the system 
and the nth environment oscillator, respectively. The total Hamiltonian is then given 
as H tot = H s + H E + Hj. In the weak coupling limit (i.e., when g <C 1), assuming an 
initially factorized state (p = Ps® Pe) and a thermal reservoir, we obtain the following 
secularly approximated master equation for the damped harmonic oscillator [U [T7j 



d , * Aft) - 7ft) , , , t , 

rfj.Ps(t) = (2a [ p s a - aa [ p s - p s aa [ ) 

A(t)+ 7 (t) , f t f v 

[2ap s a ] - a ] ap s - p s a ] a) , (2) 



where 
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A(t) =2 dt' duo J(cu) 
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(3) 



f f°° J(ui) 

7(t) = 2 / dt' dco^-^sm(ujt')sm(uj t'), (4) 
Jo Jo 2 

are the diffusion and dissipation coefficients. N(u) = {e w l kBT — 1) _1 is the average 
number of reservoir thermal excitations, with ks the Boltzmann constant and T the 
reservoir temperature, and J{oS) is the spectral density of the environment defined, in 
the continuum limit, as 

J( W ) = a 2 V-^%- Wli ), (5) 
^— ' m n u n 

with m n the masses of the environmental oscillators. No Markovian approximation was 
done in obtaining equation (J2]) so this equation describes accurately the environment 
memory. The memory effects of the reservoir are contained in the time-dependent 
coefficients, given by equations ([3]) and (jl]). It has been shown that performing the 
secular approximation does not affect considerably the system dynamics, provided that 
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Figure 1. Scaled spectral distribution I = I/(g 
r = uj c /luo = 0.1. 



■ 2 ksT) as a function of cj = uj/uiq for 



we are in high T weak coupling limit and restrict our attention to the parameter regime 
r = u c /u <C 1 [18J, where co c is a reservoir parameter introduced in Sectional Now we 
can proceed to the description of the reservoirs. 

3. Modeling the reservoirs 

The spectral densities I examine are of the form 



The exponential cutoff with cutoff frequency oo c is introduced to eliminate divergencies 
in the u> — » oo limit. The parameter s appearing in equation ([6]) is a constant that 
can acquire values < 1, 1 or > 1, corresponding to the so called sub-Ohmic, Ohmic and 
super-Ohmic spectral densities, respectively. In this paper I consider some examples and 
fix the value of s to 1/2, 1 and 3. The three cases describe different physical contexts 
(see discussion in [7]). 

The spectral distribution gives full information on the reservoir properties and it is 
given by I(co) = J{oj) \N{uj) + |] . In high T this approximates to = J(u)k B T/u. 
The parameter r = u c /uq characterizes the overlap of the system frequency with respect 
to the reservoir. In this paper I focus on a parameter region where r 1. From previous 
studies |?1 [19] we know that this is the regime where non-Markovian effects are most 
pronounced. In figure 1 all three reservoirs are shown for r = 0.1. With this choice 
of the parameters the time scale of the system dynamics, given by t$ = 1/u}q, is much 
shorter than the relevant reservoir time scales tr = l/u c . For this reason one would 
expect the non-Markovian effects to be strong. 

4. Wigner function dynamics for an initial squeezed state 




(6) 



The Wigner function is a phase space representation of the density matrix of the system 
and is thus a complete description of the state of the system. The solution of the master 
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equation d2J), based on algebraic properties of superoperators presented in [5], is used 
here because it provides a solution in terms of the quantum characteristic function 

Xt(0 = e -A rW |^ Xo[e r W /2 e -^ L (7) 

which can then be directly used to obtain the Wigner function. Here xo is the initial 
quantum characteristic function, £ is a complex variable, and 

r(t) =2 / jitjdh, (8) 



A r (t) = e- r(t) / e^Afa) dt v (9) 



We want to study initial squeezed coherent states. These states are obtained from the 
vacuum by operating first with the displacement operator D = exp(ao<^ — a^a) to obtain 
a coherent state, and then with the squeezing operator S = exp(l/2(;a:a 2 — z*a^ 2 ), 

\a ,z)=S(z)D(a o )\0). (10) 

Here z = se~ % ^ is the squeezing parameter. 

The quantum characteristic function for an initial squeezed coherent state is 

X o(£) = exp[-i|£C s - Ce^S s \ 2 + z(£*a* + £a )], (11) 

where C s = cosh(s) and S s = sinh(s). By taking the Fourier transform of the quantum 
characteristic function ([7j), with the help of equation ffTT]) . we obtain the Wigner function 
for the initially squeezed and displaced vacuum state («o = 0) with squeezing angle 
6 = 



W t (a) = Mexp 



(12) 



where 



-r(t) -2s 

(Ax) 2 (t) = A r (t) + 2 6 (13) 

-r(t) 2s 

(Aym = A r (t) + e -^— (14) 

are the variances of the dimensionless quadratures x = (a+a^) / v^2 and y = — i(a— a^) y/2 
and M is a time-dependent normalization constant [20] . 

At the initial time t = we have (Ax) 2 (0) = e~ 2s /2 = a 2 /2 and (Ay) 2 (0) = 
e 2s /2 = 1/2<T 2 . In the following I consider the squeezing dynamics in the x quadrature. 
The squeezing condition is (Ax) 2 < 0.5. For larger variances the state is not squeezed. 
An example of the squeezing dynamics for the Ohmic, super-Ohmic and sub-Ohmic 
reservoir is shown in figure 2. 

Oscillations in the squeezing for non-Markovian times appear for Ohmic and sub- 
Ohmic reservoirs. They are induced by the temporarily negative values of the coefficients 
A(t) and 7(t) which are known to occur at high T and for r 1, as in the case 
considered here. In this region, virtual exchanges of excitations between the system and 
the reservoir characterize the dynamics (THj. These virtual processes are causing the 
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Figure 2. Variance of quadrature x for a) Ohmic, b) super-Ohmic and c) sub-Ohmic 
reservoir. Markovian result is given by dashed line. In b) the Markovian dynamics 
is so fast it appears here as an almost vertical line. The parameters arc g = 0.1, 
uj c /uj a = 0.1, a 2 = 0.1 and k B T/uj Q = 1.5 x 10 3 . 

non-Markovian oscillations in the squeezing. The Markovian behaviour of the squeezing 
is plotted for comparison as a dashed line. The Markovian squeezing is obtained by 
inserting into equation (iTHjl the Markovian values of the coefficients A(t) and 7(£), 
which are given as A M = ttJ(uj q ) and 7 m = izI(u)q)/2. 

The parameters affecting the squeezing oscillations are the reservoir type, reservoir 
temperature and initial squeezing a 2 . The reservoir type is mostly responsible for the 
form of the curve, while the temperature and initial squeezing merely shift the curve 
with respect to the squeezing-non-squeezing border, (Ax) 2 = 0.5. 

The variance (Ax) 2 exhibits very similar behaviour when coupled to Ohmic and 
sub-Ohmic reservoir but is different in character for the super-Ohmic case. The super- 
Ohmic reservoir induces non-monotonic dynamics in the variance, but for the choice of 
parameters used in the figures the squeezing oscillations do not appear. From the form 
of the curves we see that for all parameters, the super-Ohmic reservoir leads to a non- 
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squeezed state in the shortest time. Oscillations in the quadrature variance persist only 
as long as the decay rate A(i) attains negative values. For the super-Ohmic reservoir 
and r < 1 the decay rates attain positive values faster than the other two reservoirs [7] . 

The differences between the squeezing for the super-Ohmic reservoir and for the 
Ohmic and sub-Ohmic reservoirs can be traced back to the reservoir properties. By 
comparing the figure 1 of the reservoirs with the squeezing plots in figure 2 we can 
easily come to the same conclusion as in [7J when considering the heating function. 
Namely, that the low frequency part of the spectrum I(oj) affects strongly the non- 
Markovian features. This means that when the spectrum is more tightly confined to the 
low frequencies (with respect to uo), longer lasting oscillations both in squeezing and 
heating are present. 

5. Conclusions 

In this paper I have studied the non-Markovian squeezing dynamics of an initially 
squeezed coherent state. Analytic expressions for the variances of the quadratures 
depend on the reservoir structure. The effect of different reservoirs was studied for 
the system coupled to differently structured reservoirs. Depending on the reservoir type 
and other parameters, the initially squeezed state showed non-Markovian oscillations 
between squeezed and non-squeezed states. 

The quantum Brownian motion model describes a large number of different physical 
systems j2U [22j [23J [2H [25] . Thus, understanding the short-time dynamics of this system 
contributes to a wide variety of experimental and theoretical scenarios. 

Reservoir engineering techniques allow in principle to tune the reservoir parameters, 
especially the parameter r, in order to reach the oscillatory regime r< 1 here illustrated 
[26J. The squeezing oscillations are connected to oscillations in the width of the 
Wigner function [16], which can be measured with homodyne detection. Therefore 
the experimental detection of the non-Markovian squeezing dynamics is, in principle, 
possible. The fact that non-Markovian dynamics are different for differently structured 
reservoirs may prove useful when considering implementing quantum devices with 
different physical setups, e.g. trapped ions or solid state materials. In any case 
understanding the effects of the environment on a quantum system is important also 
from a purely theoretical point of view. 
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